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CALCULATION PROCEDURE FOR TRANSIENT HEAT TRANSFER TO A COOLED 
PLATE IN A HEATED STREAM WHOSE TEMPERATURE 
VARIES ARBITRARILY WITH TIME 
by James Sucec* 

Lewis Research Center 

SUMMARY 

Solutions for the surface temperature and surface heat flux are found for a flat 
plate, convectively cooled from below, when the temperature of a fluid flowing over it 
in a laminar, constant property, slug- flow fashion is varying arbitrarily with time. A 
specific solution is given for the important case, from the practical standpoint, of the 
fluid temperature at the plate leading edge varying in a linear fashion with time between 
two specified temperatures. The solution to the linear case, after being properly gen- 
eralized, also serves as the basis for the rapid, approximate treatment of transients 
induced by any fluid temperature variation with time. This is done by approximating 
the actual fluid temperature variation at the plate leading edge by a series of ramps or 
steps, thus giving relatively simple functions for the surface heat flux and temperature 
responses. Several example cases are worked to illustrate the application of the known 
solutions for the step and the ramp to more general fluid- inlet temperature variations. 
Finally, a limited comparison with some experimental data is made. 


INTRODUCTION 

This work predicts the unsteady surface heat flux and temperature of a plate when 
the fluid passing over it has its temperature at the plate leading edge varying arbitrarily 
with time. As such, it represents an extension and generalization of the work reported 
previously in reference 1. 

* Associate Professor of Mechanical Engineering, University of Maine, Orono, 
Maine; NASA & ASEE Summer Faculty Fellow at the Lewis Research Center in 1972 
and 1973. 



Transients in a plate are often induced by a controlled change with time of the tem- 
perature of the fluid flowing over the plate as opposed to the case of a controlled change 
in the plate temperature or heat flux. When the fluid induces the transient, which the 
plate must then respond to, the temperature distributions in the fluid and the plate are 
mutually coupled, and one has a so-called ’’conjugate’' problem. This plate could be an 
idealization of a rocket motor wall during startup, a nuclear reactor component, a re- 
cuperative heat exchanger during startup and shutdown, or a regenerative heat ex- 
changer either in the transient unsteady domain or, the ultimate periodic unsteady time 
domain. The primary physical motivation for the work, however, is the transient in- 
duced in gas turbine blades and vanes as a result of startup, shutdown, or changes in 
the steady- state power level of an already operating engine. 

A review of much of the previous associated work in transient, forced- convection 
heat transfer is given in reference 1. The most pertinent of the works reviewed are 
references 2 and 3, in which heat transfer at a stagnation point due to a sudden change 
in free- stream temperature, with the wall temperature field mutually coupled to the 
fluid temperature field, is solved using approximate integral techniques. In reference 4 
the dynamic response of thin heat exchanger walls is treated by the Laplace transforma- 
tion under the assumption of quasi- steady conditions with a constant surface coefficient 
of heat transfer. Reference 5 treats the case of a transient in a pipe wall (infinite thick- 
ness) due to a change in the fluid inlet temperature but restricts the analysis to quasi- 
steady conditions and to the use of a constant surface coefficient of heat transfer between 
the wall and the fluid. One of the problems dealt with in reference 6 is the unsteady heat 
transfer to and temperature of a sphere when it is suddenly subjected to a fluid with a 
different temperature flowing in a Stokesian manner about the sphere. The entire 
sphere is lumped in the space coordinates, and its temperature is coupled to that of the 
fluid. A Laplace transformation of the governing equations, followed by numerical in- 
version, yields the heat- transfer response curves for the sphere. Quasi- steady solu- 
tions are also presented, and their range of validity discussed. In reference 7 a fluid 
whose temperature is varying sinusoidally with time enters a duct in a slug flow fashion, . 
that is, with a uniform velocity distribution. The duct walls are thin plates, insulated 
on their outside surfaces and interacting with the fluid on their inside surfaces. After 
neglecting, axial conduction in plate and fluid, lumping the plate temperature trans- 
versely, and using the energy balance equation on the plate as a fluid boundary condition, 
a technique somewhat like the method of complex temperature (ref. 8) is used to effect 
a solution by separation of variables. Reference 9 also deals with the case of a sinus- 
oidal fluid temperature at the entrance of a duct, but with nonparticipating walls. 

The unsteady surface heat flux and temperature of cooled gas turbine blades and 
vanes, caused by a turbine inlet temperature varying in some prescribed, though com- 
pletely general, fashion with time, is the subject of the present investigation. Using a 
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number of idealizations, a solution to the problem is found which allows study of the 
general heat transfer behavior for laminar flow over blades and vanes (as well as any of 
the. applications previously mentioned) during this type of transient. 


ANALYSIS 


Arbitrary Variation of Fluid Temperature at x .= 0 with Time 


Consider the plate of, thickness b (fig. 1) with its bottom surface subject to a cool- 
ant at constant temperature T c and with a constant surface coefficient of heat transfer 
h between the coolant and the plate bottom. (Symbols are defined in appendix A. ) The 
fluid passing over the plate has its temperature at x = 0, the plate leading edge, chang- 
ing with time t in some prescribed fashion (fig. 2). For t < 0, T = T . The response 

V 

of the surface heat flux and temperature to this arbitrary variation of fluid temperature 
at x = 0 is required. 

As in reference 1 the flow is considered to be laminar, low-speed, constant pro- 
perty, two-dimensional planar, and of the thin thermal boundary layer type. In addi- 
tion, a steady slug flow velocity profile is used, plate and fluid axial conduction is 
neglected, and the plate temperature is to depend on only x and t-(i.e. , the plate’s 
temperature is lumped in the y direction). In attempting to find the response to an 
arbitrary variation with time of the fluid temperature at x = 0, one first seeks the solu- 
tion for a step change in the fluid temperature which then, by virtue of the linearity of 
the governing equations and boundary conditions, can be used with Duhamel's theorem 
to yield the more general result. The mathematical statement of the problem for a step 
change in fluid temperature from T c to T Q is 


3T 9T 3 2 T 

— + u — = a 

3t 3x -2 


9y 


at t = 0 x > 0, y>0T=T c at x = 0 t > 0, y > 0 T = Tq as y — t > 0, x>0 
T remains finite. The exact analytical solution to this equation, subject to the side 
conditions shown, was found in reference 1 using double Laplace transformations. The 
solution is now presented in terms of the more convenient variable (for use with 


Duhamel's theorem) a = T - T 


with cr 0 = T 0 -T c : 


2 \ 
a = ct 0 u(t - 1) ^erf (Y) + e 2T?Y+T? {erf [e(r - 1) + tj + Y] - erf (77 + Y)}j 


( 1 ) 
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where the unit step function is 


u(t - 1) =< 


ro 

i+i 


for t < 1 
for r — 1 


( 2 ) 


The appearance of the unit step function is the mathematical manifestation of the 
physics of the plane wave phenomenon, namely, that at position x the plate cannot re- 
spond to what has happened to the fluid at x = 0 and time = t until the fluid, that was 
at x = 0 at time = t, reaches the position x. The nondimensional variables appearing 
in equation (1) are defined in reference 1 and repeated here for convenience. 




(3) 


Using Duhamel's theorem (ref. 8) along with equation (1) when ct q = 1 and /3 = <r gives 
the transient temperature distribution within the moving fluid for arbitrary variation of 
the fluid temperature at the plate leading edge with time as 


a 



/3(t?,€, Y, r - X) 


dg x = Q( X > 

dX 


dX 


(4) 


where X is a dummy variable for the nondimensional time t. 

Since the plate surface temperature is usually of primary interest, Y is set equal 
to zero and the function /3 is inserted to yield the result for the plate surface tempera- 
ture when the fluid temperature variation is arbitrary, namely, 



{erf [e(<p - X) + tj] - erf (rj)} 


d W X > 

dX 


dX 


(5) 


where y = t - 1. The treatment of the unit step function during the integration is shown 
in appendix B. Defining a surface heat flux parameter F as 
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and noting that 



F = - 



allows one to use Duhamel's theorem to get the analogue of equation (5) for the surface 
heat flux. Thus, the surface heat flux when the fluid temperature at x = 0 is varying 
arbitrarily with time is given as 



[e 2 (<p-X) 2 +277€(<p-X)] 


+ 




(6) 


When using either general expression (eq. (5) or (6)) one must keep straight the fact that 
X is the dummy variable for nondimensional time r, not for cp . Both integrals, of 
course, are to be interpreted in the Stieltjes sense. 

As a quick check of equation (5) and illustration of its use, it will be applied to the 
case of a step change in temperature of the fluid at t = 0 from T £ to Tq and then 
compared with the previously derived result (eq. (15) of ref. 1). At an abrupt change in 
CT x=0 a * ^0’ i nte £ ra l equation (5) becomes equal to the bracketed portion of the 
integrand evaluated at Xq times the change in ct x _q that occurs at Xq. For the step 
change at t = 0 one has Xq = 0 and Aa x _Q = <7 x _q - ct”_q, where the plus refers to the 
value just after t = 0 and the minus signs refers to the value just before t = 0, hence 
Ao x _q = Oq. Equation (5) thus becomes 

cr 2 

— = e 71 {erf(e<p + rj) - erf(?j)} for q> > 0 (7) 

CT 0 


Rearrangement shows that this equation is identical to equation (15) of reference 1 as it 
should be. 
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Solutions for Generalized Step and Ramp in Fluid Temperature at x = 0 


In general, for an arbitrary T X=Q as a function of time, the integrals of equations 
(5) and (6) often cannot be found analytically. But T X=Q can be approximated by a se- 
quence of steps and ramps to any desired accuracy, and equations (5) and (6) do admit 
analytical, relatively simple solution functions for both the step change and the ramp 
change in fluid temperature at x = 0. Hence, response functions for any arbitrary fluid 
temperature at x = 0 can be constructed, approximately, by the proper combination of 
the response functions for the generalized step and generalized ramp as will be ex- 
plained later. 

For the response to a generalized step in fluid temperature where T x _ Q changes 
from T x _q to T x=q at time tj, which corresponds to r. = X, equations (5) and (6) 
yield the solutions (from ref. 10) 


A(T w i = AT i eV ( erf " V + *l] " erf fa)} 


( 8 ) 


for cp > and 


*?- 


■AT^ 


- t 2 (^-'r i ) 2 +277e(^-'r i )| r 2 , 

; J + \ir tie 71 (erf [e(<p - 7.) + 77 ] - erf (i))j 


(9) 


for (p > r., where ATj = T x=0 - T~ =0 the subscript i on the Act w and F refers to 
the time at which the step occurs, and the superscript s indicates that these are re- 
sponses to a step as opposed to a ramp. 

The generalized ramp, in which T x=Q changes in a linear fashion from T. at time 
to Tj at time tj, is depicted in the portion of figure 2 bounded by the dashed refer- 
ence lines. The response functions for this case are important not only for their role 
when the fluid temperature varies arbitrarily, but also because often during a startup, 
shutdown, acceleration, or deceleration of a gas turbine engine, the turbine-inlet tem- 
perature, as a first approximation, can be represented as varying linearly with time. 
The equation of the generalized ramp is 


£=o - T i + < T j - T i ) (rri) * 4 * 


( 10 ) 


Using this in equations (5) and (6) leads to the following response functions for the ramp 
(as presented in ref. 10). 
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( 11 ) 



Tj) erfc (tj) + i 1 erfc £e(<p - t.) + 7 ?] - i 1 erfc 



for 


< (p < T. 


Aa ^' — rj~ { € ^ T i ~ erfc ^ + iJ erfc [ € ^ 


Tf) + 


i* erfc 




for 



( 12 ) 




- r\e{<p - Tj) erfc ( 77 ) 


+ Tjji 1 erfc ( 77 ) - i 1 erfc |e(<p - t.) + Tj]jj (13) 


for Tj < cp < Tj and 



for cp > Tj. The subscript i denotes the time at which the ramp begins, the super- 
script r is the response due to the ramp, the prime on the r indicates the validity of 
the expression in the (p domain given by r. < (p < Ty and i 1 erfc ( 77 ) is the first re- 
peated integral of the error function and is defined and tabulated in reference 11. (A 
representative portion of the details of the integrations of eqs. ( 5 ) and ( 6 ) subject to 
eq. (10) is given in appendix C. ) 


Solutions Obtained by Combining Steps and Ramps 

Equations ( 8 ), (9), and (11) to (14), when appropriately combined, give the transient 
surface temperature and heat flux due to a prescribed fluid temperature variation at 
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x = 0, when it is approximated by a series of steps and/or ramps. As an illustration, 
consider the fluid temperature variation at x = 0 shown as a solid line in figure 3. 

Also shown in figure 3 are three dashed lines that form a step and two ramps which ap- 
proximate fairly closely the actual fluid temperatures as a function of time. (No dashed 
line is shown beyond time t g because it coincides with the solid line there. ) Hence, the 
flux at time t, say, when tj < t < t 2 is given by the following sum: 


F = (Eq. (9) with Tj = 0 and AT. = Tj - T c ) + (Eq. (14) with Tj = 0, r- = Tj, and Tj 


- T i = T 2 - Tj) + (Eq. (13) with Tj = Tj, t. = r 2 , and T. - T. = T Q - T g ) 

The details of the procedure for the case depicted in figure 3 are illustrated in ap- 
pendix D. 

RESULTS AND DISCUSSION 


Transient Due To Change From One Steady State Operating Level to Another 


In the analysis section and appendix D, the procedure for handling the more general 
cases of fluid temperature variation with time are discussed, and the example shown in 
figure 3 is worked in some detail. Now consider the case depicted in figure 4, where 
T x=0 = Tj for t < 0 and then linearly increases to T Q at time t = tj and is held at 
Tq thereafter. 

The response of the surface temperature and flux consists of the proper combining 
of equations (8), (9), and (11) to (14). The contribution to a w and F for t < 0 can be 
considered to have resulted from a step change in fluid temperature from T c to Tj a 
long time ago. Hence, this portion of the response is given by equations (8) and (9) as 
(p (interpreted here as the time elapsed since the step function changed T x Q from T 
to Tj) goes to infinity. Thus, the contributions to the wall temperature and flux for 
t > 0, but due to the history of the fluid before t = 0, are 


2 

Aor w = ( T i " V eT? erfc to) 


(15) 


and 


F = 



erfc (tj) 


(16) 


For 0 < t < tj the right hand side of equation (11) would be added to equation (15) and 
the right hand side of equation (13) would be added to equation (16) with r. set to zero 
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and t. = Tj with T. - T. = T Q - Tj. These sums give the responses in the time do- 
main between 0 and tj. For t > tj the right hand sides of equations (12) and (14) 

would be added to equations (15) and (16), respectively, to yield surface temperature and 
heat flux for values of time t > tj. The solution is now complete. Division of the final 
equations for the a w and F by T Q - T £ to yield (T w - T C )/(T Q - T c ) and Q w , re- 
spectively, gives, as the solution functions for the surface temperature and surface heat 

flux, for the case shown in figure 4, the following: 




erfc ( €<f + tj) 


- i* erfc (rj)j (18) 


for 0 s < Tj 


T„. - T„ /t, - T \ 2 


JM erfc (i,) J T ° ' TAe " f - - 1 


T 0 - T c \ T 0 - T c/ 


^ T 0 - T c 


\ * 

i\£L{, 

c/ €T l 1 


eTj erfc (rj) + i erfc {ecp + 77 ) 


i 1 erfc je(<p - Tj) + rjjj (19) 


for <p ^ Tj 


erfc w * ( £ [erf ^ + 


‘ T 0 ‘ T c/ 


r j) - erf tj] 


+ tj€ (p erfc (tj) + Tjji 1 erfc {ecp + tj) - i 1 erfc (tj)]| (20) 


for 0 ^ <p < Tj, and 


% ^ ^ erfc M + (2 H {e<t> + ” > " erf ■ Ti> 

+ TJ j + TjeTj erfc (t?) + Tjji 1 erfc [e<p + tj] - i 1 erfc |e(<p - Tj) + tj]|^ (21) 


fOr (P S- Tj. 
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Equations (18) to (21) are plotted in figures 5 and 6 for a value of exj = 2. 00 and 
for (Tj - T C )/(T 0 - T c ) = (T Q - Tj)/(T 0 - T c ) = 0. 50. The exj value is (as noted in 
ref. 10) a measure of the steepness of the ramp that takes T x _ Q from Tj at t = 0 to 
T 0 at t = tj. The curves in figures 5 and 6 are plotted for two values of tj: tj = 0 cor- 
responding to an insulated plate (not cooled) and tj = 1 corresponding to a rather well 
cooled plate. The plate, of course, does not respond to the change in the fluid temper- 
ature at t = 0 until t = 1, hence, the reason for plotting the time parameter as 
e(r - 1). It is seen from figure 5 that, for both cases, tj = 0 and t? = 1, the wall tem- 
perature increases monotonically from the original steady-state temperature to the final 
steady- state temperature corresponding to the new fluid temperature at x = 0, T Q . In 
figure 6 the surface heat flux continually increases until the displaced time (r - 1) at 
which the fluid temperature at x = 0 ceases its variation with time. Beyond this time 
the surface heat flux decays to its final steady state value. 


Comparison With Data 

A comparison of the predictions of the theory presented in this work with some ex- 
perimental data would be desirable. However, only an extremely rough comparison 
could be made because of incomplete data and, in particular, data from a case that did 
not satisfy the major restrictions of the idealized model. The measured data was from 
an air cooled vane (described in ref. 12) during an acceleration in a laminar flow region 
of the vane surface. Because of a nonslug velocity profile, axial and spanwise conduc- 
tion, transverse temperature gradients, and x-dependent h Q in the experiment, one 
would not expect even the initial or final steady- state nondimensional temperatures to be 
predicted well. Hence, an attempt was made to over-ride some of these effects by com- 
paring (T w - T c> i n itiap/( T w, initial " T c, initial since this P arara eter is unity at the 
start of the acceleration for both the measured data and the theoretical solution. The 

subscript ’'initial” pertains to values of T c and T w just before the acceleration. 
During the transient in the turbine-inlet temperature, which was approximated by a 
single ramp for the analysis, one also encounters time varying T , h , and u 
These parameters are assumed constant and are, in addition to the aforementioned 
effects, therefore not taken into proper account in the analysis. Also present in the 
experiment was extreme property variation and a variable thickness wall. To make a 
rough comparison, arithmetic averages of T £ , u w , and the fluid properties were used 
along with the value of h £ at the end of the acceleration. Considering the vast differ- 
ences between the experimental conditions and the constraints on the analytical model, 
the agreement between theory and experiment is at least encouraging. In figure 7 the 
measured wall temperature response curve is plotted against time at x/L = 0. 40 for 
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an acceleration in which the turbine-inlet temperature changes from 922 to 1644 K 
(1200° to 2500° F). Also plotted is the analytical curve, which comes from equations 
(18) and (19). The qualitative behavior of the nondimensional temperature excess ratio 
just referred to was predicted by the theory. The percentage error between the two was 
about 26 percent at low times and about 12 percent at the final steady state. The mani- 
fold differences between the theoretical model and the experiment preclude any conclu- 
sive explanation of these percent differences. However, a rough calculation indicates 
that the neglect of transverse temperature differences (the lumped in the y direction 
assumption) is an important difference between the theory and experiment. 


SUMMARY OF RESULTS 

Expressions for the surface temperature and surface heat flux are found for a plate, 
convectively cooled from below, while the temperature of a fluid passing over the top of 
the plate is changed arbitrarily with time at the plate leading edge. 

1. Simple, easy to use functions were presented that allow a solution to the problem 
of fluid temperature varying arbitrarily with time by approximating the arbitrary fluid 
temperature variation by a sequence of ramps or steps. 

2. The mechanics of properly combining these functions to yield plate temperature 
and heat flux is demonstrated by means of examples. 

3. Equations are given for the important case of a ramp in fluid temperature at 

x = 0 causing T x _q to change from temperature Tj at time t = 0 to temperature Tq 
at t = t^ and held at Tq thereafter. 

4. A limited comparison of the predictions of the idealized model with experimental 
results was made for a case that did not satisfy many of the restrictions of the model. 
This yielded a rough qualitative similarity of the predictions and the data. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 4, 1975, 

505-04. 
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APPENDIX A 


SYMBOLS 


b 

C. 


f ; 

r* 

F i 

F S 

k c 

i 1 erfc (z) 

k f 

% 

% 

T 

T c 

T x=0 

t 

u(r - 1) 

U oo 

V 

X 

Y 

y 

a f 

P 

€ 

P 


plate thickness 

constant- pressure specific heat 
surface heat flux parameter, 


i w 


V 7rQ?£S 


surface heat flux response function defined by eq. (14) 
surface heat flux response function defined by eq. (13) 
surface heat flux response function defined by eq. (9) 
coolant- side surface coefficient of heat transfer 

✓■CO 

first repeated integral of error function, / erfc (z') dz’ 
thermal conductivity of fluid flowing over upper surface of plate 
instantaneous nondimensional surface heat flux 
instantaneous surface heat flux 
temperature 
coolant temperature 

instantaneous fluid temperature at x = 0 
time 

unit step function, equal to 0 for r < 1, and to +1 for r^l 

free -stream velocity 

variable defined by eq. (C5) 

space coordinate along plate 

nondimensional y coordinate defined in eqs. (3) 

space coordinate perpendicular to plate 

thermal diffusivity of fluid flowing over plate 

value of a when ctq = 1 

nondimensional variable defined in eqs. (3) 

density 
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X dummy variable for time 

a temperature difference, T - T c 

act^ temperature response function defined by eq. (12) 

ao^ temperature response function defined by eq. (11) 

A<7 w temperature response function defined by eq. (8) 

i 

77 nondimensional variable defined in eqs. (3) 

r nondimensional time defined in eqs. (3) 

Tj nondimensional duration of ramp that begins at t = 0 

cp shifted nondimensional time, r - 1 

Subscripts : 

f properties of fluid flowing over plate 

i index 

j index 

w wall 

0 ultimate steady- state fluid condition at x = 0 

°° free- stream conditions 

Superscripts : 

r response to ramp after it has ended 

r' response to ramp before it has ended 

s response to step in fluid temperature at x = 0 



APPENDIX B 


DEVELOPMENT OF EQUATION (5) 

To find a w for an arbitrary variation of fluid temperature at x = 0, against time, 
when the wall is initially at T , Y is set equal to 0 in equation (4), and /3 is given by 
the right side of equation (1) with Y = 0, a Q = +1, and t replaced by r - X. Doing this 
yields 




/■ 

*/0 


(r - X - 1) jerf [e(r - X - 1) + ?j] - erf (rj)| 


d( W X) 

dX 


dX 


(Bl) 


However, u(t - X - 1) = +1 only for r - X - 1 > 0 by definition of the unit step function; 
hence, X < t - 1. Thus, since X is the integration variable and values of X > r - 1 
lead to a zero integral, this is taken into account by making the upper limit equal to 
r - 1, which has been defined as (p. This change obviates the need for explicit appear- 
ance of the unit step function in equation (Bl) and thus it becomes 


2 ff , 

% -•’’■/ {' 

Jn 


i d<T v-n( x ) 

erf [e(<p - X) + rj] - erf (rj)i — — dX 

J dX 


which is equation (5). 


(B2) 
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APPENDIX C 


DEVELOPMENT OF EQUATION ( 11 ) 


Inserting equation ( 10 ) into equation ( 5 ) when 
ment, 


T. 


i < * < T i 


gives, after rearrange- 




[e(<P - X) + rjJ - erf 



Define the integral as I: 


I = 



\e{<p - X) + 77] - erf 



dX 


(Cl) 


(C 2 ) 


The second term of equation (C 2 ) can be immediately integrated to give the following ex- 
pression for I: 


Define 



erf [< z(<p - X) + 77] dX - (<p - r.) erf (77) 


(C 3 ) 


Ij =/ erf [e(<p - X) + 77] dX 


(C4) 


To facilitate the integration in equation (C 4 ), the following variable change is made: 


v = e(<p - X) + 77 


(C 5 ) 


Therefore, 


dX = - 


dv 


e 


and this variable change causes the limits on the integral to be mapped as follows: r. 
corresponds (from sq. (C 5 )) to e(<p - Tj) + 77 = v and cp corresponds to v =77, so 
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equation (C4) becomes 


i rv i MV' 7 i) +r ? 

I. = - - I erf (v) d v = - / erf (v) dv (C6) 

e •'e(<p-Tj)+Tj e *tj 

Adding and subtracting unity to the integral of equation (C6), noting that 1 - erf (v) = 
erfc (v), and integrating the constant term gives 


i M<P- tJ+tj 

1^ = <p - t. — J erfc (v) dv (C7) 

€ *TJ 


Rearranging gives 




r°° 

erfc (v) dv - I erfc (v) dv 

•^((p-T.J+rj 


(C8) 


However, the remaining integrals are simply the first repeated integral of the error 
function (ref. 11), and thus equation (C8) becomes 


Ij = cp - T i + - ji 1 erfc [e(<? " T j) + tj] - i 1 erfc (rj)j 


(C9) 


Substituting equation (C9) into (C3) and then equation (C3) into (Cl) yields the final re- 
sult 




je(<p - 7-p erfc (rj) + i 1 erfc | ~e(<p 


r 0 + v] 


i 1 erfc 



(CIO) 


for Tj < <f < Tj. 


This is equation (11). 
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APPENDIX D 


MECHANICS OF SURFACE TEMPERATURE FUNCTION SOLUTION FOR 
FLUID TEMPERATURE VARIATION SHOWN IN FIGURE 3 

The solid line of figure 3 shows the variation with time of the fluid temperature at 
x = 0. The surface temperature response as a function of time is desired. Rather than 
using the actual fluid temperature against time in equation (5) and then numerically in- 
tegrating the result, the approximate procedure stressed in this report will be used. 
Hence, the actual fluid temperature variation is approximated by a step and two ramps 
as shown by the dashed lines in figure 3. Now equations (8), (11), and (12) must be 
properly combined to yield the surface temperature as a function of time. 

First consider the time period 0 < t < tj, which corresponds to 0 < t < Tj and, 
because cp = r - 1, also to 0 ^ <p < Tj. In this time period the fluid temperature at 
x = 0 is still on the first ramp; hence, the plate has only felt the effect of the step 
change and the portion of the first ramp that has occurred. Hence, its temperature re- 
sponse is the sum of equations (8) and (11), properly interpreted. (Note that eq. (12) is 
not yet needed because <p is less than Tj not greater. Also, no response function for 
the second ramp is needed because the second ramp has not yet occurred. ) Since the 
step occurs at t = 0, it follows that r. = 0 for the step and that ATj = Tj - T c in 
equation (8). Also, the first ramp begins at t = 0 or Tj = 0 and ends at tj, which 
gives Tj = u^tj/x and (Tj - T.) = (T 2 - Tj). Utilizing this information in equations (8) 
and (11) and then adding them gives the wall temperature response for 0 ^ <p < Tj as 

2 /t - T \ 2 

<7 W = (Tj - T c )e 1? [erf (e<p + -q) - erf (?j )] + f— -j e 71 [tcp erfc (rj) 

+ i 1 erfc (ecp + 77) - i 1 erfc ( 7 ])] (Dl) 

In the time interval t^ — t < t 2 > corresponding to Tj < <p < T 2 (where Tp = u co t 2 /x) 
the surface temperature response consists of the response to the step, which is simply 
the first term of equation (Dl), the response due to ramp 1, and the response due to the 
portion of ramp 2 that has occurred. Since ramp 1 is over, cp > Tj, which means 
<p > Tj for ramp 1; and equation (12) applies with r. =0, Tj = t j, and T. - T. = 

Tg - Tj. For the response to ramp 2, equation (11) applies, since for this ramp 

T i = T l» T j = t 2> and T 1 ~ V < t 2‘ ^ 1so for ram P 2 > Tj - Tj = Tq - T 2 . With these 
substitutions, the sum of the first term of equation (Dl) (eq. (8) really) and equa- 
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tions (12) and (11) give the surface temperature in the nondimensional time domain 
Tj ^ <p < T 2 as follows: 


2 Xt - T \ 2 

a w = ( T 1 " T c> e7? [ erf ( €< P + V) - erf (t?)] +( — -J ZeTj erfc ( 77 ) + i 1 erfc (e<p + tj) 

eT l / 1 


i 1 erfc \e{(p - Tj) + tj] j + — je(<? - Tj) erfc (tj) 


.1 

+ 1 


erfc | \{(p - Tj) + rjJ - i 1 erfc (tj)| 


(D2) 


For times t > t 2 corresponding to <p > T 2 , the surface temperature response is 
due to the step at r = 0, the first term of equation (D2), the ramp between 0 and Tj, 
which is now complete, the second term of equation (D2), and the ramp between Tj and 
t 2 , which is also complete. This last contribution is given by equation (12) with 

T j " T i = T 0 ” T 2’ T i = T V and T j = T 2’ Nothin S else has to be added beyond t 2 since 
there the fluid temperature at x = 0 is a constant. Hence, adding the first two terms 
of equation (D2) to equation (12) gives the surface temperature in the time domain 
<p > Tj as 


2 /p T \ 2 

a w = (Ti - T c )e T? [erf (e<p + tj) - erf (tj)] +f— ^ -j e 77 jeTj erfc (tj) + i 1 erfc (ey + tj) 


- i* erfc \^{<p - Tj) + 


•)} 


e(r 2 - Tj) erfc (tj) 


.i 1 


erfc 


\z(<p ~ Tj) + Tj] - i 1 erfc \e{<p - t 2 ) + tj]| 


(D3) 


The procedure for finding the surface heat flux function F is essentially the same 
as that which led to equations (Dl) to (D3), except that equations (9), (13), and (14) will 
be used instead of equations (8), (11), and (12). To demonstrate the steady-state form 
of a w for this case, and also as a partial check on equation (D3), let (p — «> in equa- 
tion (D3). Since erf (°°) = +1, 1 - erf (tj) = erfc (r\), and i 1 erfc (°°) - 0, equation (D3) 
reduces to, after cancelling and combining terms, 
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2 

<x w - ( T o ‘ T c ) e7? erfc W as (p - «o 

This is also what equation (7) reduces to as <p — «>, and this is as it should be, 
since at long times after the last disturbance is complete (ramp 2 in this case as shown 
in fig. 3), the wall "forgets" the nature of the disturbances which changed T x _q from 
T c to Tq. Hence the steady- state solution should be, and is, independent of the way 
in which T x _q goes from T c to Tq. 

It should now be apparent how the procedure used here on a step and two ramps is 
to be used when any number of steps or ramps are chosen to approximate an arbitrary 
fluid temperature at x = 0 as a function of time t. 
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Figure L - Wall flux against shifted time for a ramp change of inlet temperature. Ramp 
steepness, - 2.00; nondimensional temperature ratios, (Tx - T c )/(Tq - T c ) - 

tf(r T i )/(T (r T c ) “ a5a 



Figure 1. - Comparison of measured and predicted wall temperature variation with time at x/L • a 4a 
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